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Wc present a rneasurcnicnt of the average multiplicity of mesons in S", _B" and B^ meson 
decays. Using 17.6 fb^^ of data taken at the T(45) resonance by the BABAR, detector at the PEP-II 
e+e" storage ring at the Stanford Linear Accelerator Center, we reconstruct <j) mesons in the K^K~ 
decay mode and measure B{B — > 0X) = (3.41 ± 0.06 ± 0.12)%. This is significantly more precise 
than any previous measurement. 

PACS numbers: 13.25.Hw 
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I. INTRODUCTION 

The large data sample collected by the BABAR de- 
tector provides an excellent opportunity for a signif- 
icant improvement to the existing measurements of 
the average (/) multiplicity in B meson decay at the 
T(45') resonance. This quantity, which is convention- 
ally denoted B(B 4>X), was previously measured 
by CLEO as (2.3±0.6±0.5)% [H and by ARGUS as 
(3.90±0.30±0.35)% % These two measurements dis- 
agree at the 1.8(7 level, leading to a large error on the 
Particle Data Group average (3.5 ± 0.7)% 0. The 
OPAL collaboration has measured the average multiplic- 
ity B{b (j)X) = (2.82±0.13±0.19)% at the Z° pole. 
This latter measurement is sensitive to 6-hadron decays 
that are not accessible at T{AS) experiments, including 
h baryons and, in particular, the meson. 

An improved measurement of B{B 4iX) can lead to 
improved measurements of the oscillation frequency. 
The primary decay modes of the B^ meson contain 
mesons, which often (18% 0) produce a meson in their 
decays. Due to this high rate, B^ decays into mesons 
are a prime decay chain for 5° oscillation searches. An 
important input to such searches is the knowledge of 
the background arising from non-strange B meson de- 
cays into (j) mesons. 

Given the large size of the BABAR data sample, this 
measurement is limited by systematic errors. As a re- 
sult, this analysis is designed to minimize these system- 
atic errors. Minimal selection criteria are applied, and 
efhciencies and backgrounds are evaluated directly from 
data where possible. The measurement is performed in 
(/)-momentum intervals to minimize the systematic effects 
that may be introduced by differences between the </> mo- 
mentum spectrum in data and simulation. 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 

The data used in this analysis were collected by the 
BABAR detector at the PEP-II storage ring. We use 
17.6 fb^^ of data taken at the T{AS) resonance (on- 
resonance) and 4.1 fb^^ of data taken at a center-of-mass 
energy 20MeV below the BB threshold (off-resonance). 
The latter sample is used for the subtraction of the non- 
BB component (continuum) in the on-resonance data. 
These data samples were taken between January and 
May 2002. Additional data, consisting of 3.5 fb~^ of on- 
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resonance data and 1 fb~ of off-resonance data taken 
under different running conditions, are used for verifica- 
tion of the result. 

A detailed description of the BABAR detector is pre- 
sented in Ref. . The components of the detector most 
relevant to this analysis are described here. Charged- 
particle tracks are reconstructed with a five-layer, double- 
sided silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH) with a helium-based gas mixture, placed 
in a 1.5-T solenoidal field produced by a superconduct- 
ing magnet. The resolution on pT, the charged-track 
momentum transverse to the beam direction, is approxi- 
mately {Spt/pt? = (0.0013 (GeV/c)-VT)^ + (0.0045)2. 
Charged particles are identified from the ionization en- 
ergy loss (AE/dx) measured in the DCH and SVT, 
and the Cherenkov radiation detected in a ring-imaging 
Cherenkov device. The efficiency for identifying true 
kaons exceeds 80% over most of the momentum range of 
interest, while the probability for a pion to be misidenti- 
fied as a kaon is less than 3%. 

We use Monte Carlo samples of 7(45") ^ B°B^ and 
B'^B~ decays, corresponding to twice the expected num- 
ber of B mesons in the data sample, to study our selection 
efficiency. The B-meson decays are simulated accord- 
ing to previously measured branching fractions which 
account for approximately 60% of all B decays. The 
remaining 40% are modeled by JETSET JJ, while pre- 
venting any enhancement of the first 60%. The detector 
response in these samples is simulated with the GEANT4 
program and cross-checked with control samples in the 
data. 



III. EVENT AND CANDIDATE SELECTION 

Events are selected if at least three tracks are found 
and the measured total energy is at least 4.5 GeV. In 
order to suppress continuum background, events are re- 
jected if the ratio of the second-to-zeroth order Fox Wol- 
fram moments |9] (-R2) is higher than 0.25. This require- 
ment rejects 62% of the off-resonance data, while retain- 
ing 78% of the simulated B (j>X events. 

The selection of — > K^K" candidates requires 
two oppositely-charged tracks that satisfy 0.1 < pr < 
10 GeV/c, have at least 12 hits in the DCH, are consistent 
with originating from the primary interaction point, and 
satisfy kaon identification criteria based on dE/dx mea- 
surements and Cherenkov radiation. Tracks are assigned 
a kaon mass hypothesis and neutral two-track combina- 
tions are formed. Candidates are selected if their invari- 
ant mass is in the range 1.004 < ttikk < 1.036 GeV/c^. 
This mass window is equivalent to about 4.5 standard 
deviations on either side of the nominal </> mass, where 
the RMS spread in the ttikk distribution is due to both 
the natural </> width and the detector resolution. This 
relatively large acceptance is chosen to reduce the effect 
of a possible mass resolution difference between data and 
Monte Carlo, at the expense of signal-to-background sig- 
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nificance. 

A total of 471, 941 (34, 900) candidates survive these 
selection criteria in the on- (off-)resonance sample. 



IV. BACKGROUND ESTIMATION 

Two sources of background to the B — > (pX signal are 
considered: random combinations of tracks that pass the 
selection (combinatorial background) and true (f) mesons 
that do not originate from BB events (continuum back- 
ground). Because the reconstruction efficiency depends 
on the momentum of the 0, these backgrounds are sub- 
tracted separately in 16 bins of (f) momentum. 

We first remove the continuum background from our 
signal by subtracting the mxK distribution obtained in 
the off-resonance sample from that in the on-resonance 
sample, scaled by the ratio of the luminosities of the two 
samples. This scale factor is calculated by comparing 
the number of e'^e~ — » //"""/i" events in the two sam- 
ples. The center-of-mass momenta of (j) candidates in 
the off-resonance data are scaled by the ratio of on/off- 
resonance beam energies to account for the slightly dif- 
ferent momentum spectrum of the continuum component 
in the on-resonance sample. This procedure explicitly ac- 
counts for all backgrounds from physics processes other 
than T{4S) production as their cross sections are almost 
identical at the two energies; it also accounts for beam- 
related backgrounds, as the running conditions were very 
similar. 

We next subtract the combinatorial background to ex- 
tract the number of (j) mesons. This background is es- 
timated by fitting the mass distribution in sideband re- 
gions well away from both the signal and the KK thresh- 
old. The ranges 0.989 < tukk < 1-002 GeV/c^ and 
1.04 < rriKK < 1.1 GeV/c^ were chosen and the function 
{rriKK — 2mi^)" ■ {b + c- tiikk + d ■ m]^^ + e ■ m]^^) was 
used. This function provides a good description of the 
phase space in the vicinity of a threshold; a fit to the 
combinatoric ttikk spectrum (removing true 4> mesons) 
in Monte Carlo gives = 104 for 95 degrees of freedom. 

Figure ^ shows the ttikk distributions of the on- 
resonance and luminosity-scaled off-resonance data sam- 
ples for all (j) momentum bins combined; the fitted com- 
binatorial background shapes are overlaid. The on- (off-) 
resonance fit has — 86 (75) for 61 degrees of free- 
dom. With estimates of the combinatorial background 
shape, the signal is extracted by subtracting the back- 
ground. The resulting signals are shown in the lower 
part of Figure n In the on-resonance sample, we observe 
(2.349 ±0.007) X 10^ candidates, and the corresponding 
number in the off-resonance sample is (1.95 ±0.02) x lO'*. 
These plots and numbers are representative only; they 
are not used in the signal extraction. 

Erroneous estimation of the signal in data and simu- 
lation could arise from a number of sources. The phase- 
space function may be unable to describe the background 
shape correctly (see Monte Carlo fit above). There may 




m„„ (GeV/c') 



FIG. 1: Top: Invariant mass distributions of candidates pass- 
ing all selection requirements except that for the mass. The 
solid histogram shows candidates in the on-resonance data 
sample, while the dashed histogram shows the off-resonance 
sample, scaled to the luminosity of the on-resonance data. 
The fitted combinatorial background is overlaid on both his- 
tograms as a dotted curve. Bottom: Resulting signal after 
combinatorial background subtraction. Again, the solid his- 
togram represents on-resonance candidates, and the dashed 
histogram represents the luminosity-scaled off-resonance can- 
didates. Note that the vertical scale is displaced from zero. 
The bump visible in the lowest bins of these plots is a thresh- 
old effect, which is well understood, and the fit range was 
chosen so that this does not contribute to the measurement. 



be (j) mesons with a reconstructed mass in the sideband 
regions, which would affect the background estimation 
and therefore the signal yield. The Monte Carlo may 
not correctly model the signal lineshape (in particular, 
the fraction of 4> mesons outside the signal mass region), 
leading to a systematic uncertainty on the efficiency. We 
consider each of these sources below. 

We vary the fitting procedure in a number of ways in 
order to test the robustness of our background estima- 
tion. We replace the third order polynomial in the above 
function with both a second order polynomial in ttikk 
and an exponential term exp(& ■ itikk)- These changes 
result in a 0.15% and 0.65% change in the number of cj) 
candidates respectively. 

To account for the possibility that the reconstructed 
(f> mass extends into the sideband regions, we vary the 
upper bound of the region excluded from the fit, raising 
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it from 1.04 to 1.06 GeV/c^ (while keeping the signal 
region as defined above). The largest difference in the 
number of (f) mesons in the signal region in the above 
variations is 2.4%. 

Finally, we look at the fraction of candidates, af- 
ter background subtraction, outside the signal mass re- 
gion. We count the number of candidates in the range 
1.036 < rriKK < 1-05 GeV/c^ and calculate the ratio of 
this number to the number of candidates in the signal 
region. This ratio is found to be 2.4% in data and 2.6% 
in Monte Carlo. This yields a difference of 0.2% on the 
number of (j) candidates between data and Monte Carlo. 

We take the largest difference in all the above tests 
(2.4%) to be the systematic uncertainty associated with 
the combinatorial background subtraction and signal se- 
lection. 



V. SELECTION EFFICIENCY 




^ ^ 

1 1.02 1.04 1.06 1.08 1.1 

"^KK (GeV/c ) 

FIG. 2: Mass distribution of the K+ K- , K+T- and K-T+ 
samples. All selection criteria were applied in the K'^K~ 
sample but one track is not required to pass the kaon selection 
in the K^T^ samples. 



The kaon identification efficiency is extracted from 
data to avoid the systematic errors associated with Monte 
Carlo-based determinations. To do this, candidates are 
constructed from two-track combinations with at least 
one track passing the kaon identification criteria. This 
is done for positive and negative tracks separately to ac- 
count for a possible asymmetry. Three subsamples of the 
data are defined: K^K~ , where both tracks have passed 
the kaon identification requirements, and K^T^ , where 
only one track is required to pass the kaon selection. The 
same off-resonance subtraction and R2 requirement is 
made in defining these samples as for the standard se- 
lection. The kaon identification efficiency is then given 
by the ratio of the number of 4> mesons reconstructed in 
the K~^K~ sample to the number in the K^T^ samples: 
eif± = N(j)KK /N(f>x^j'±. Figure 13 shows the invariant 
mass distribution of the K^K^ , K^T^ and K^T^ sam- 
ples over the entire momentum region. 

Studies of data and Monte Carlo samples show that 
the kaon identification efficiency is not constant through- 
out the momentum range of our <j) sample, but can be 
described by one constant efficiency values below = 
1.2 GeV/c and another above this value. This step in effi- 
ciency is caused by the transition from dB/dx-based par- 
ticle identification at low momenta to Cherenkov-angle- 
based particle identification at higher momenta. Since 
our analysis was performed in (f) momentum bins, such 
behavior may introduce a bias in the result if the kaon 
selection efficiency is taken as a constant value over the 
entire range. We therefore extract the kaon selection ef- 
ficiency separately above and below this momentum. We 
measure £k+ = (98.6 ± 1.2)% and Ek- = (98.7 ± 1.2)% 
for p0 < 1.2 GeV/c, and eK+ = (79.0 ± 1.8)% and 
Ek- = (79.7 ± 1.9)% for > 1.2 GeV/c. 

The remaining efficiency to be estimated is £2t, that of 
finding two charged tracks that originate from a (j) meson, 
satisfy the R2 requirement, and have an invariant mass 
in the signal region (with no kaon identification require- 



ment). This efficiency is estimated from Monte Carlo. 
Since the efficiency to reconstruct a (f) depends on the 
(j) momentum, differences between the momentum spec- 
trum in data and the generated spectrum in the Monte 
Carlo sample must be considered. Therefore, the analy- 
sis was carried out separately in 16 bins of 4> momentum. 
The bins are chosen to have equal (with the exception of 
the lowest momentum range) numbers of reconstructed (f) 
mesons in the Monte Carlo. Figure|21shows the efficiency 
e2T as a function of the (f) momentum. 



uT 0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



J , , I I I 







0.5 



1.5 2 2.5 

(GeV/c) 



FIG. 3: Efficiency (e2T) as a function of <p momentum. The 
highest bin includes all entries above 1.9 GeV/c. Errors shown 
are statistical only. 
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VI. RESULTS 



The average multiplicity is calculated with the formula 



B{B (l)X) 



1 



2NbbB( 



K+K- 



16 

■E 



N' 



B.i 



(1) 



where Ng ^ is the number of (p mesons in momentum bin 
i found in the data and assumed to come from B mesons. 
This number is obtained by performing the background 
fit to the on-resonance data samples after subtracting the 
off-resonance data samples, scaled to the on-resonance lu- 
minosity. The efficiency e = e2T£K+^K- is the product 
of the reconstruction efficiency and the kaon identifica- 
tion efficiencies for each track. The quantity Nbb is the 
number of BE events in the data sample, which is mea- 
sured to be Nbb = (18.7 ±0.2) x 10^ using a technique 
described elsewhere 

Since the analysis was performed in momentum bins, 
the efficiency in each bin has very little dependence on 
the modeling of the </> spectrum, except for the lowest- 
momentum bin, which includes the tracking detection 
limit. We therefore sum the yield in the highest 15 bins 
and extrapolate the result based on the simulated spec- 
trum, so that the sum in Equation is replaced by 



16 

■ 1 



E 

1=2 



Nt 



MC,i 



y^l6 jy. 



(2) 



MC,i 



Here, iV^c % ^^'^ number of 4> mesons in the Monte 
Carlo sample in momentum bin i. 

Using S(0 ^ K+K-) ^ 0.492 ± 0.006 0, we obtain 
B{B <pX) = (3.41 ± 0.06)%, where the error is statis- 
tical only. Figure ^ shows the measured and simulated 
(j) momentum spectra in the T{AS) center-of-mass frame. 
We note that the Monte Carlo sample predicts the ob- 
served (j) momentum spectrum reasonably well. 



^ 0.04 
O 0.03 



& 0.02 
§ 0.01 
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FIG. 4: Center-of-mass 4> momentum spectra of the measured 
and Monte Carlo generated samples. The measured spectrum 
shows statistical and systematic errors combined, and the gen- 
erated spectrum is normalized to the measured multiplicity. 
The lowest momentum bin has a negative central value due to 
the large continuum component but is consistent with positive 
values due to the large error. 



TABLE L Relative systematic uncertainties. 



Source 


AS/B(%) 


Combinatorial BG fitting 


2.4 


On/Off scale factor 


0.1 


5(0 -> K+K-) 


1.2 


Nbb 


1.1 


£2T (total) 


1.9 




0.5 


Monte Carlo statistics 


0.3 


Monte Carlo 4> modeling 


0.7 


Tracking efficiency 


1.6 


Total 


3.4 



VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are associated with all of the 
variables in Equation 1. Tabled lists the various system- 
atic uncertainties, which are described in detail below. 

Two sources contribute to the uncertainty on iV^^, the 
number of (j) mesons from BB events. One is the fitting 
procedure error which is described in Section Hvl and was 
taken to be 2.4%. The other is the error on the scaling 
factor (on/off- resonance) which contributes 0.1% relative 
uncertainty on the measured average multiplicity. 

The uncertainty on ^(0 — >• K^K~) introduces a rela- 
tive error of 1.2%, while that on Nbb contributes 1.1% 
relative uncertainty on the average multiplicity. 

As £2T is obtained from Monte Carlo, differences be- 
tween data and Monte Carlo give rise to systematic un- 
certainties. Several components contribute to this sys- 
tematic uncertainty: 



• There is a systematic uncertainty related to the ex- 
trapolation of the 15-bin yield to the full result, 
due to our limited knowledge of the (j) spectrum in 
B decays. Since only about 60% of B meson de- 
cays are well understood (see Sectional, mismod- 
eling of the remaining 40% can affect the result. 
To account for this model dependence, we study 
two Monte Carlo subsamples representing extreme 
cases to make up the entire remaining 40%. The 
first subsample contains a B meson undergoing a 
two-body decay to a charm meson, with the charm 
meson undergoing a two-body decay to a (j). In the 
second subsample, the B meson undergoes a multi- 
body (greater than two) decay and the subsequent 
charm meson undergoes a multi-body decay into a 
final state that contains a (j) meson. These two cases 
yield very different kinematic distributions for the 
meson. We measure the fraction of candidates in 
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the lowest bin in each sample. We take the largest 
difference between these samples and the primary 
result systematic uncertainty. It is found to 
be 0.7%. The effect of (f) polarization was similarly 
studied and found to have a negligible impact on 
the result. 

• To establish the contribution to the systematic un- 
certainty from the simulation of i?2, the number of 
4> mesons from BB events is estimated again with- 
out the i?2 requirement, and the same procedure 
is applied to the Monte Carlo. The fraction of (f) 
mesons from B decays with R2 < 0.25 in data is 
(78.18±0.8G)% while this fraction in Monte Carlo 
is (78.00±0.09)%, in agreement within statistical 
errors. We also investigate the two decay models 
mentioned above for their effect on the R2 selec- 
tion. We find that the largest difference between 
the models and our Monte Carlo distribution is 
0.5%, and we take this difference as a systematic 
uncertainty. 

An additional test is performed by examining dif- 
ferent continuum-suppression variables. We study 
the angle between the </> direction in the center-of- 
mass frame and thrust axis of the event, where the 
thrust was calculated both including and exclud- 
ing the (j) candidate. These two variables are each 
used in place of R2 in order to suppress continuum 
events. We place appropriate criteria on these vari- 
ables to maintain similar efficiency to that of R2 in 
our analysis. Wo then measure the efficiency of 
these requirements in data and Monte Carlo. The 
ratio of efficiencies between data and Monte Carlo 
is found to be 0.982 ± 0.015 for the first variable 
(with the 4>) and 1.007 ± 0.015 for the other. 

• Tracking performance is studied using control sam- 
ples in data, and the track-finding efficiency is 

foimd to be accurate to within 0.8% per track. We 
therefore assign a 1.6% systematic uncertainty due 
to tracking efficiency. 

• Finally, the statistical uncertainty on e2T con- 
tributes a 0.3% systematic uncertainty. 

The use of one single kaon selection efficiency for all 
(f) momenta was compared to the use of separate values 
above and below = 1.2GeV/c. The observed differ- 
ence in the average multiplicity was 0.9%. This is below 
the statistical error on the kaon identification efficiencies, 
hence no additional error was assigned to this source. 
The statistical error on the kaon selection efficiencies is 
treated as part of the statistical error in this analysis as 
it is obtained from the same data set as our signal and 
scales appropriately. 

The above sources of systematic uncertainty are added 



in quadrature and yield a relative uncertainty on the av- 
erage multiplicity of 3.4%. 

This analysis is repeated by replacing the 16 bins 
with 6 bins of 6*0, the polar angle of the ((> candidate with 
respect to the beam axis. As the total number of events 
is exactly the same as in the analysis described above, 
this is not an independent measurement, and can only 
serve to validate the fitting procedure. The combinatorial 
background in these bins is significantly different in shape 
to that used in the primary analysis. The total yield 
of (j) mesons from B decays is found to differ by 0.88% 
from the yield in the primary analysis — well within the 
assigned uncertainty for this source. 

We also repeat the analysis using a different data set. 
We use a smaller data set from the year 2000 in which the 
detector was operating under different conditions. This 
analysis yields B{B (j)X) = (3.34 ± 0.07)% where the 
error is statistical only, entirely consistent with our pri- 
mary result. 



VIII. CONCLUSION 

By selecting two identified oppositely-charged kaons 

from a sample of T{AS) data and subtracting the combi- 
natorial and continuum background, we measure the av- 
erage multiplicity of 4> mesons in B meson decays. Our 
measurement of B{B (bX) = (3.41 ± 0.06 ± 0.12)% is 
consistent with both previous measurements at the 1.5ct 
level, although it is significantly more precise. 
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